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Abstract 

We compute up-down asymmetries for atmospheric neutrino induced electron and muon 
events in the SuperKamiokande detector for the following neutrino oscillation models: (A) 
maximally mixed — v T or — u s , (B) maximally mixed — v e , (C) threefold maximal 
mixing between v e , and v T , and (D) neutrino oscillations via Equivalence Principle or 
Lorentz invariance violation. We emphasise the role of different momentum cuts. 



The solar |TJ] and atmospheric neutrino observations and the LSND || exper- 

iment provide strong evidence that neutrinos have nonzero masses and oscillate. Many 
specific neutrino mass scenarios have been proposed as explanations for some or all of 
these results. These explanations will be put to stringent tests in the next few years as 
more data are collected. 

In this paper we will focus on ways to experimentally test several different proposed 
explanations of the atmospheric neutrino anomaly. One good way to discriminate between 
these possible explanations is through the up-down asymmetry, with respect to zenith 



angle, of the detected charged leptons in experiments such as SuperKamiokande |T0|||TT[] . 



This idea was introduced in Ref. [01 and studied in detail in Ref . [ITOfl in the context of 



the Exact Parity Model solution to the atmospheric neutrino anomaly \12\. The up- 



down asymmetry idea was independently considered by the authors of Ref. [ITT] in the 



case of six possible neutrino oscillation models. As emphasised in Ref. [11], the sign and 



energy dependence of the up-down asymmetry will provide a useful diagnostic tool in 
differentiating the various explanations of the atmospheric neutrino anomaly. Reference 



TTJ focussed on the energy dependence of the up-down asymmetry. A complementary 
approach is to study the asymmetry with various cuts in momentum [|irj[|13j. This is 
useful because the Kamiokande and SuperKamiokande collaborations have divided their 
event sample into two classes (sub-GeV and multi-GeV) on the basis of momentum cuts. 
In the future, improved statistics may lead to the event sample being further subdivided 
with respect to momentum. The purpose of this paper is to study the up-down asymmetry 
expected at SuperKamiokande with various cuts in momentum for a variety of neutrino 
oscillation models. 

It is a priori interesting to consider various cuts in zenith angle as well as energy when 
analysing the up-down asymmetries of charged lepton events in SuperKamiokande [|l(J • A 
class of measures of the up-down asymmetry can be defined by introducing the quantities 
yl and y^, where 

Here N~ v {N^ v ) is the number of charged electrons produced in the detector with zenith 
angle cos 6 < — r\ (cos 6 > rf). N^ v are the analogous quantities for charged muons. 
The zenith angle G is defined so that G = corresponds to downward travelling and 
G = 7r to upward travelling charged leptons. Different choices for rj correspond to different 
cuts in zenith angle. (Note that Ref.[pH| considered the rj = case in detail for several 
different models, while Ref. JlTJ studied different choices for rj for the Exact Parity Model.) 
Since y^ are ratios of fluxes, they will be approximately free of the systematic errors 
arising from uncertain absolute fluxes (±30%) as well as uncertain fi/e flux ratios (± 
a few percent). For symmetry reasons, ~ 1 is expected in the absence of neutrino 
oscillations. If neutrino oscillations are present, however, then y^ ^ 1 is expected. 
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In the calculations to be presented later, we actually use the related quantities Y£ , 
where 



yri — 



yV 


osc 




no~osc 



(2) 



The numerators are the predictions for in the models, while the denominators are the 
same quantities in the absence of oscillations. We will focus on the case 



7] = 0.2 



(3) 



because our experience with various values for rj in Ref . |T(J suggests that this choice leads 
to the largest effects without compromising too much in regard to statistics. Using rj = 0.2 
also allows comparison of our results with the preliminary data from SuperKamiokande. 
The cases we consider are 

(A) Maximal 



m 



or v^ — Vg mixing [O. These two scenarios are indistinguishable 
as far as up-down asymmetries are concerned provided that matter effects can be neglected 
(which is approximately the case for the region of 5m 2 considered in this paper). The 
~ v s scenario is, in particular, motivated by the Exact Parity Model in the region 
of parameter space where v e — v' s oscillations can be neglected for atmospheric neutrinos 
(Ref.fjHJ focussed on the parameter space region where v e oscillations are also important). 



(B) Maximal — v e mixing. The Acker-Pakvasa three-flavour model 0] is essentially 
indistinguishable from this scenario for atmospheric neutrinos. 

(C) Threefold maximal mixing [pj]][[il|] amongst z/ e , and u T in the region of parameter 



space considered in Ref. |TS|. 

(D) Massless neutrinos with violation of the Equivalence Principle or breakdown of Lorentz 
invariance The special case of exactly maximal u e — oscillations 0] is considered 
for defmiteness (sin 2 29 can actually be as low as 0.8 phenomenologically) . 



The Cardall-Fuller three-flavour scheme |15] is numerically similar to case A with the 
parameter choice 5m 2 ~ 0.3 eV 2 . For such a large 5m 2 there are no expected up-down 
asymmetries for the sub-GeV or multi-GeV SuperKamiokande data. 

Consider sub-GeV neutrinos first. In the water- Cerenkov Kamiokande and Super- 
Kamiokande experiments, the sub-GeV neutrinos are detected via the charged leptons £ a 
(£ a = e or /j,) produced primarily from the quasi-elastic neutrino scattering off nucleons in 
the water molecules: u a N — > £ a N' (a = e,fi). The event rate is calculated by integrating 
the product of the differential neutrino flux, the scattering differential cross-section, the 
energy efficiency function for the detector and the relevant neutrino oscillation probability 



with respect to energy and angular variables (see, for example, Ref.|10[ for more details). 
In our numerical work we have used the differential cross section in Ref.[21|. The differen- 
tial flux of atmospheric neutrinos without geomagnetic effects is given in [22|, but we have 
used the differential flux which includes geomagnetic effects f23|. (For other atmospheric 
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neutrino flux calculations, see Ref . |[24|| . ) The detector energy efficiency function can be 
found in Ref. . 

In principle, the charged lepton event rates for multi-GeV neutrinos require a more in- 
volved computation, because several scattering channels must be considered. Fortunately, 



as noted in Ref.flll]], the up-down asymmetries are generally insensitive to these details. 
In practice, all one really needs to calculate up-down asymmetries reliably is a reasonable 
estimate for a function which is effectively a product of the cross-section and the energy 
efficiency function. We have estimated this function from the SuperKamiokande Monte- 
Carlo calculations of event rates in the absence of oscillations. We numerically checked 
the dependence of the asymmetries on variations in this "effective cross-section". Cases 
A, B and C vary by only a few percent, while case D turns out to be uncertain to about 
the 10% level. 

In Figures 1-6 we plot Y e °„ 2 for each of the cases A, B, C and D, for three different 



momentum cuts, as a function of the relevant parameter |26| . For case A the parameter is 



the difference of squared masses between and either v T or v s . For case B it is the 5m 2 
between z/ M and v e . In case C it is the larger of the two independent 5m 2 's, while for case D 
it is the usual measure of either Equivalence Principle or Lorentz invariance violation. In 
Figs. 1 and 2 we have utilised the standard Kamiokande sub-GeV momentum cuts, 0.2 < 
p^/GeV < 1.5 and 0.1 < p e /GeV < 1.33, for muons and electrons respectively. Figures 3 
and 4 show the corresponding plots when the alternative cuts, 0.5 < p^/GeV < 1.5 and 
0.5 < p e /GeV < 1.33, are employed. Figures 5 and 6 pertain to multi-GeV events. We 
now discuss some of the important qualitative features of the graphs: 

Most of these curves show that the asymmetry depends weakly on the neutrino oscilla- 
tion parameter (5m 2 or Sv/2) over a significant range (almost an order of magnitude). This 
can be understood as follows: Neutrinos with cos8 > 0.2 travel distances 15 ~ L/km ~ 
75, whereas neutrinos with cos 8 < -0.2 travel distances 2600 ~ L/km ~ 12, 700. Thus 
there will be a large range of 5m 2 or Sv/2 where the oscillation length is such that neu- 
trinos with cos 8 > 0.2 do not have time to oscillate, whereas neutrinos with cos 8 < —0.2 
oscillate and are averaged. This weak dependence is interesting because it effectively pro- 
vides a prediction within each model for the asymmetries which is reasonably insensitive 
to parameter choice and hence avoids the potential problem of fine-tuning. (Note that 
the plateau feature in the asymmetries is also there when r\ = is used instead of rj = 0.2. 
The point is that, due to the geometry of the situation, the neutrino flux from the inter- 
mediate regime —0.2 < cos© < 0.2 forms a small enough fraction of the total flux that 
the qualitative plateau phenomenon persists.) 

A comparison of Figs. 1 and 3 shows that case A can be more clearly distinguished from 
the other cases by adopting the alternative sub-GeV momentum cut. Case A is in turn 
clearly different from the no-oscillation Y — 1 case. Cases B and D can be distinguished 
from case C through electron asymmetries, particularly using the alternative sub-GeV 
cut, according to Figs. 2 and 4. Although case C is not clearly different from the no- 



3 



oscillation situation for electron asymmetries, one can return to the muon asymmetries 
of Figs. 1 and 2 to obtain a clear differentiation. The most problematic differentiation for 
the sub-GeV sample is evidently that between cases B and D. If 5m 2 or Sv/2 is sufficently 
large [>~ 1CT 2 eV 2 or (km GeV) -1 ], then this discrimination can be more easily made. 
Fortunately, the multi-GeV cases displayed in Figs. 5 and 6 show a clear difference 
between cases B and D unless 5m 2 or 5v/2 is quite small. 

Preliminary atmospheric neutrino results from SuperKamiokande can be found in 
Ref.[]27|]. The asymmetries Y^ 2 can be estimated from the graphs presented for the sub- 



GeV case with the usual momentum cut and for the multi-GeV case. The results are, 

Y°- 2 (sub - GeV) ~ 1.19 ± 0.10 (stat.) Y^ 2 (sub - GeV) ~ 0.80 ± 0.07 (stat.) (4) 
and 

Y°- 2 (multi-GeV) ~ 0.80±0.15 (stat.) Y°- 2 (multi - GeV) ~ 0.46±0.07 (stat). (5) 

These numbers can be compared to Figures 1, 2, 5 and 6. These figures show that 
while case A fits the preliminary data best, the still significant statistical errors preclude 
definitive conclusions. The Cardall-Fuller model fits the data least well. 

In summary, up-down asymmetries in the charged leptons induced by atmospheric 
neutrinos in the SuperKamiokande detector are important quantities that, with improved 
statistics, will be capable of clearly distinguishing the various proposed solutions to the 
atmospheric neutrino anomaly. 
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Figure Captions 

Figure 1. Y^ 2 with momentum cut 0.2 < p^/GeV < 1.5 for case A (solid line), B (dashed- 
dotted line), C (dashed line) and D (dotted line). Note x = 5m 2 (eV 2 ) for cases A, B 
and C while x = 5v/2 (km~ 1 GeV~ 1 ) for case D. Also shown (straight dotted lines) is the 
preliminary superKamiokande result (Eq.(|J) together with 1 sigma statistical errors. 

Figure 2. Y^ 2 with momentum cut 0.1 < p e /GeV < 1.33 for the cases A (solid line), B 
(dashed-dotted line), C (dashed line) and D (dotted line). Also shown (straight dotted 
lines) is the preliminary superKamiokande result (Eq.([|) together with 1 sigma statistical 
errors. 

Figure 3. Same as Figure 1, except that the momentum cut 0.5 < p^/GeV < 1.5 has 
been taken. 

Figure 4. Same as Figure 2, except that the momentum cut 0.5 < p e /GeV < 1.33 has 
been taken. 

Figure 5. Y^ 2 for the SuperKamiokande multi-GeV sample. Also shown (straight dotted 
lines) is the preliminary superKamiokande result (Eq.(^) together with 1 sigma statistical 
errors. 

Figure 6. Y" e 2 for the SuperKamiokande multi-GeV sample. Also shown (straight dotted 
lines) is the preliminary superKamiokande result (Eq.(^]) together with 1 sigma statistical 
errors. 
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